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INITIAL RESULTS FROM THE VOYAGERS 1, 2 MAGNETIC FIELD EXPERIMENTS 


M. H. Acuna, L. F. Burlaga, 

R. P. Lapping and N. F. Ness 
Laboratory for Extraterrestrial physics 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 


1. Introduction 


The aim of this paper is to survey some of the initial results 
obtained with the Goddard Space Flight Center magnetic field 
experiments aboard voyagers 1 and 2. We consider the period from 
launch in August 1977 to March 1978, when the spacecraft moved from 1 
AU to 2.5 AU. A plot of the trajectories is given in the paper of 
Belcher et al., which appears in these proceedings. 

The interplanetary magnetic field observations on each spacecraft 
are made with a dual magnetometer system, described by Behannon et al. 
(1978). The Principal Investigator is N. F. Ness, vector magnetic 
field measurements are made at a rate of 16.66 samples/s (every 60 
milliseconds) , allowing a detailed study of small scale phenomena such 
as shocks. At high bit rates, all of the measurements are 
transmitted, while at the lowest data rates 1.92 second averages of 
the field are transmitted. Since the spacecraft are always in the 
solar wind, mesoscale and macroscale phenomena can also be 
investigated, the only interruptions being due to limits on the 
telemetry coverage. The complex task of data reduction is 

carried out jointly for the plasma science and the magnetic field 
experiments at GSFC, so plasma and magnetic field data are 
simultaneously available for analysis in a common format. Dr. H. 
Bridge, the Principal Investigator of the plasma experiment, and two 
of his Co-Investigators at MIT, J. Belcher and J. Sullivan, provided 
the plasma data that appear below. 

The plan of this paper is to: discuss some large-scale 
characteristics of the interplanetary magnetic field in the period 
August 1977-March 1978; describe a characteristic of a post-shock 
flow; and analyze the detailed structures of three types of 
interplanetary shocks. 

2 . Large-Scale Structure 

Hourly averages of the magnetic field intensity (B) and direction 
(x,5) are shown in Figures 1 and 2 for two consecutive solar 
rotations. The angles \ and fi are measured in heliographic (RTN) 
spacecraft centered coordinates; 6 is the latitude measured with 
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respect to the RT 
plane, (fi points 
away from the sun? T 
lies in a plane 
parallel to the 
solar equator and 
points in the 
direction of the 
sun's rotation); x 
is the longitude, 
which is zero when 
the field points 
along R and which 
increases 
counterclockwise 
when viewed from the 
north . 
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The magnetic 
field variations in 
Figures 1 and 2 are 
very different from s 

the stable, 
recurrent two sector 
pattern observed by 
Pioneers 10 and 11 
when they moved 
through the same 
heliocentric radial 
distance range in 

1972/1973 (Smith and Rms 

7 FEB 

Wolfe, 1977). Most 
significantly, 
large, recurrent 
sectors do not 
appear in the 

Voyager data. Small regions of uniform polarity can be seen in Figure 
1, but they do not recur on the following solar rotation (Figure 2) 
with the same general characteristics. There are several intervals in 
which the magnetic field does not lie near the nominal spiral 

direction (x ^ 135 or 315° a:t 1 AU for negative or positive 
polarization) . The magnetic field intensity profile in Figures 1 and 
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2 is likewise irregular and non-recurrent. (Note the logarithmic 
scale.) For example, exceptionally low intensities were observed on 
January 29 (Figure 1) , but they did not recur on the following 
rotation. Similarly, exceptionally low intensities were observed on 
February 14 (Figure 2) , but they were not observed 27 days earlier. 

In general, the magnetic field intensity varies more on SR 1976 than 
on SR 1975. 

The results in Figures 1 and 2 are representative of the 
large-scale fields observed by Voyagers 1 and 2, i.e., the magnetic 
field intensity and direction are irregular and they usually do not 
show simple, recurrent patterns. The reason that the Voyager results 
are so different from the Pioneer results is probably associated with 
the phase of the solar cycle. The Voyager observations were made 
during an increase of solar activity, when there were many active 
regions at mid-latitudes and few large, stable coronal holes near the 
equator. The corresponding Pioneer observations were made in the 
declining phase of the solar cycle, when activity was low and large, 
and stable coronal holes were present near the equator. 

3. A Post-Shock 
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Flow 

A shock and 
post-shock flow were 
observed near 1 AU 
by Voyagers 1 and 2 
and by Earth 
orbiting IMPs 7 and 
8 in the period 
September 19 to 21, 

1977, when the 
Voyager 1-sun-earth 
angle was 4.2° and 
the Voyager 
2-sun-Voyager 1 
angle was 4.3° as 
shown in Figure 3. 

The multiple 

spacecraft configuration was ideal for studying azimuthal gradients, 
i.e., in the direction normal to the earth-sun line. Behind the shock 
at Voyager 2 (0507 U.T., September 20), indicated by the abrupt 
increase in speed V and magnetic field intensity B in Figure 3, the 
speed and field intensity decreased monotonically to the preshock 
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values 20 hrs later. In contrast, behind the shock observed at • 
Voyager 1 (2130 U.T. + 50 m, September 19) the speed and field 
increased for several hours before decreasing to the preshock values. 
In other words, the post-shock speed and magnetic field intensity 
profiles at Voyager 2 resembled a blast wave or P-type shock 
(Hundhausen, 1972) , while the profiles behind the shock at Voyager 1, 
just 4.2° away, resembled a driven wave or R-type shock (Hundhausen, 
1972) . The speed profile at IMPS 7 and 8 was similar to that at 
voyager 1. 

Our interpretation of these observations is that the shock was 
driven by a piston with a thin lateral boundary. The piston passed 
Voyager 1 and IMP, but it did not extend far enough in the azimuthal 
direction to intercept voyager 2. Thus, the width of the boundary was 
less than 4.3°. 

4 . Shock Structure 

Voyagers 1 and 2 observed more than 20 shocks in the 6 months 
interval under consideration. Here, only three shocks will be 
discussed: a perpendicular forward shock (B perpendicular to the shock 

A/ 

normal, and the shock moving away from the sun relative to the 
unshocked plasma); an oblique forward shock; and a quasi-parallel 
reverse shock (moving toward the sun relative to the unshocked 
plasma) . Table 1 gives the basic parameters just before and after 
each shock. 

Perpendicular 
forward shock . 

Figure 4 shows the 
shock mentioned 
above, observed at 
Voyager 2 at 0507 UT 
on September 20, 

1977 at 1.05 AU. 

The data are 1.92s 
averages obtained in 
the low bit rate 
mode . Two 
characteristics of 
the data in this 
figure identify the 
event as a 
quasi-perpendicular 
shock; 1) the large, 
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abrupt increase in the magnetic field intensity such that = 

Table 1 and Burlaga, 1971), and 2) the absence of any 
significant change in the direction of Given v, , V_ , and B /B = 

•L ^ C, A. 

'^2'^^1 mass conservation equation to obtain the shock 

speed, viz., U = 400 km/s, or U-V^^ == 60 km/s relative to the solar 
wind ahead of the shock. 

A full least squares fit to a subset of the Rankine-Hugoniot 

conservation relations (Lepping and Argentiero, 1971) yielded the 

normal orientation as = 29.1°, 6^ = -11.2° with a rest frame shock 

speed of 382 km/sec. The corresponding angle between the best fit 

pre-shock field and the normal was 88° thus confirming the 

perpendicular character of this shock. The "quality factor" normally 

associated with this calculation was 0.74 (where unity corresponds to 

an excellent fit) indicating a good fit. Using a time delay 

calculation between the sudden commencement time (1142 U.T. on 19 

September) at earth and the shock passage at Voyager 2 reasonably 

consistent results were obtained: v_ = 310 km/s, x =27°, 6 = -2°. 

s n n 

The slightly lower shock speed (19%) may be explained by local shock 
curvature in the Earth-Voyager region, consistent with the limited 
exten-t of the associated driver plasma. 

The structure of the shock is simple. There are no fluctuations 
upstream; the intensity increased in a very short time in the main 
shock transition, and behind the shock there are random fluctuations 
of very small amplitude. 

Oblique forward 
shock . An oblique 
forward fast shock 
was observed at 
Voyager 2 at 1655:10 
UT on February 2, 

1978, High 
resolution magnetic 
field data (60 ms 
samples) in Figure 5 
show an abrupt 
increase in B, which 
is characteristic of 
a shock. 

Simultaneous 
increases in the 
plasma parameters 
were also observed 
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(see Table 1) , confirming that the event was probably a forward shock. 
The direction of the magnetic field changed slightly across the shock, 
indicating that the shock was oblique. One can use the coplanarity 
theorem to calculate the shock normal (e.g., see Hundhausen, 1972). 

The angle between the shock normal and the upstream magnetic field was 
47°. The speed of the shock with respect to the upstream solar wind 
was U-V^ = 102 km/s. 

The structure of the oblique shock shown in Figure 5 is somewhat 
more complicated than that of the perpendicular shock described above, 
but it is not turbulent. Wave packets with quasi-sinusoidal , high 
frequency fluctuations appear upstream of the shock transition, and 
lower frequency fluctuations appear downstream. 

Figure 6 shows 
power spectra in the 
,i.ange (0-8) Hz for 
the undisturbed 
upstream region (A) , 
the region just 
ahead of the shock 

(B) , and a region 
behind the shock 

(C) . Spectrum C 
shows appreciable 
low frequency power 
with respect to the 
background (Spectrum 
A) , but no spectral 
peak is evident. 

Spectrum B shows a 
distinct peak at ^ 

1.2 Hz in the 5-second interval ahead of the shock, corresponding to 
the quasi-periodic -fluctuations evident in Figure 5. 

A polarization analysis for the quasi-periodic fluctuations just 
ahead of the shock showed that the fluctuations are nearly circular RH 
polarized waves propagating 40° with respect to the 

ambient magnetic field. This is shown in the holograph of Figure 7 
which gives a plot of the perturbation vector in the maximum variance 
plane; the wave vector, Jj,, and the component of B parallel to are 
normal to the plane of this figure and point into the paper. Waves 
such as these have been observed upstream of the earth's bow shock. 
Greenstadt et al . (1975) and Fairfield and Feldman (1975) identified 
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them as whistler 
waves. Taking 
advantage of the 
successive inward 
and outward motion 
of the bow shock, 
they showed that the 
bow shock whistlers 
were standing waves. 

Their technique 
cannot be applied to 
an interplanetary 
shock, which locally 
moves in just one 
direction. On the 
other hand, the 
fast, uniform, known 
motion of 
interplanetary 
shocks offers the 
advantage that scale 
lengths can be 
determined more 
precisely. 

Quasi-parallel , 
reverse shock . 

Figure 8 shows an 
unusual event , 
recorded by Voyager 
2 at 0920 UT on 
January 29, 1978. 

There is a large 
abrupt decrease in 
density and proton 
thermal speed 
(temperature) and a 

slight decrease in magnetic field intensity coincident with an abrupt 
increase in bulk speed. This is the signature of a reverse, fast 
shock (e.g., see Burlaga, 1971). 

High resolution (1.92s) magnetic field data in Figure 9 show a 
small change in magnetic field direction (X.) and a small increase in 
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intensity (B) at j; 

0919 UT, which is 
the main transition 
in the shock. The 
small size of ax and 
AB indicates that 
the shock is 
quasi-parallel (for 
a parallel shock AX 

= AB = 0 ) n 

The structure 
of this shock 
differs appreciably 
from that of the 
perpendicular and 
oblique shocks 
discussed above. 

Upstream (late 

times ) , there are low frequency, irregular fluctuations extending > 
500,000 km from the shock. A preliminary power spectral and 
polarization analysis reveals these fluctuations to be left hand, 
elliptically polarized waves with characteristic periods of 20-25 
seconds. The main shock transition is turbulent, with large amplitude 
fluctuations in both intensity and direction; these fluctuations 
extend to high frequencies, as indicated by the high variance of the 
60 ms measurements of the magnetic field components in Figure 9. 
Downstream (early times--in the shocked plasma) there are high 
frequency, turbulent fluctuations extending > 250,000 km behind the 
shock, power spectral and polarization analyses indicate that these 
fluctuations are right-hand, nearly circular-polarized waves with 
characteristic periods of 6 to 12 seconds, the period decreasing as we 
approach the main shock transition. 

5 . Summary 

Some early results of the Voyager 1 and 2 magnetic field 
measurements between 1 AU and 2.5 AU were surveyed. 

In the macroscale data (> 27 days) it was noted that the magnetic 
field intensities and sector patterns were complex and variable, in 
contrast to the simple and stable patterns oberved in the same region 
(1^0 to 2.5 A.U.) by Pioneer 10 and 11 during 1974. 

On a smaller scale, the flow configuration behind the September 
14-21, 1978 shock was investigated using magnetic field and plasma 
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data from Voyagers 1 and 2 and from IMP 7 and 8. The shock was driven 
by a piston and the azimuthal boundary of the piston was found to be 
very thin at 1 A.U. 

The structures of three different types of interplanetary shocks 
were analyzed. A perpendicular forward shock was found to be thin and 
laminar, with almost no fluctuations upstream and downstream. An 
oblique forward shock was found to be thin, with whistler wave packets 
upstream and higher intensity turbulent fluctuations downstream. A 
quasi-parallel, reverse shock was found to be turbulent, with 
irregular, LH , elliptical-polarized fluctuations extending >500,000 
km upstream, and with the nearly circular-polarized waves extending > 
250,000 km downstream. 


TABLE 1 

Shock Characteristics 


Event at Voyager 2 


®2 

"l 

"2 

^1 

^2 

^1 

^2 

Sept. 20, 1977 

3.1 

6.8 

9.1 

19.8 

2.4 

4.5 

341 

381 

Feb. 2, 1978 

2.3 

4.3 

1.7 

4.0 

1.4 

3.8 

335 

390 

Jan. 29, 1978 

0.7 

1.2 

0.5 

1.2 

1.0 

5.4 

370 

330 


(B is in gammas, n in cm T in 10‘^°K, and v in km/s) 
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ALFVEN WAVES AND ALFVENIC FLUCTUATIONS IN THE SOLAR WIND 


K.W. Behannon and L.F. Burlaga 
Laboratory for Extraterrestrial Physics 
Godda rd ^Space Flight Center 
Greenbelt, Maryland 20771, USA 


Abstract Recent measurements have confirmed that fluctuations in the 
bul k velocity of the solar wind are correlated with fluctuations in the 
magnetic field more than one-third of the time and that these observed 
"Alfvenic" fluctuations have several characteristics in common with 
A1 f ven waves. They are not "pure" Alfven waves because there are fluc- 
tuations in magnetic intensity, and at times tangential discontinuities 
are also present. Because the fluctuations in ^ are three-dimensional, 
they are most consistent with being general AlTven waves, although they 
are usually found to be "nearly" transverse. The assumption that the 
observed waves are plane waves may not be valid over distances > 0.01 AU. 
A better theoretical understanding is needed with respect to mode-cou- 
pling, stability, and propagation. 

1. INTRODUCTION 

Although the i nterpl anetary magnetic field has a spiral configuration 
on average, there are always fluctuations in it on scales of an hour 
or more. In general, these fluctuations are complex and are the result 
of many different physical processes. There is a relatively simple sub- 
set of these fluctuations which occurs more than one third of the time 
and is characterized by a high correlation between fluctuations in the 
bulk velocity ^ and the magnetic field B. These are called Alfversic 
f 1 jctuati ons . Alfvenic fluctuations are important for several reasons: 
they scatter cosmic rays; they modify heat transport; they influence 
geomagnetic activity; they may carry information on the fundamental 
acceleration process of the solar wind; and their relative simplicity 
offers the hope that their physical nature can be understood, thus con- 
tributing to a better understanding of a fundamental plasma. 

The subject of Alfven waves has been discussed in several reviews 
during the past few years (Hollweg, 1974, 1975, 1978; Barnes, 1974; 

Volk, 1975 ) which were written from a theoretical point of view. A 'com- 
prehensive, balanced review of the theory and observations of Alfven 
waves and Alfvenic fluctuations was recently written by Barnes (1978). 

In this present review we attempt to give a concise account of the basic 
types of Alfven waves that are possible and a summary of the observations 
of Alfvenic fluctuations that have been made using both magnetic field 
and plasma data. We shall tend to stress the ways in which Alfvenic 
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fluctuations (AF's) differ from pure transverse Alfven waves (TAW's), 
in the hope that it will stimulate fresh thinking and new progress in 
solar wind and wave studies. We shall also show the ways in which AF's 
resemble TAW's. 


2. BASIC TYPES OF ALFVEN WAVES 


"Alfven wave" is a general term covering several different possible 
types of magnetic field and plasma fluctuations. Necessary conditions 
for an Alfven wave of any type include the following. 

1. B(t) can be written as the sum of an unpeirturbed component B and 
the wave perturbation B,: 

B = B^ + B, . (1) 

Note that <B> = B„ if and only if <B.> = 0, and <B.> is not neces- 
sarily zero in general (here < > denotes the average). 

2. The wave moves along B„ at the speed 

•^0 


where p is the density. 

3. There is a perturbation in bulk velocity, V,, which is related to 
B, by the equation 

Ii = ± (3) 

A 

The minus sign refers to waves moving in the B^ = B^/B^ direction. 

4. The waves satisfy a 1 i near wave equation. Both large amplitude 
Alfven waves (Bj<^B^) and snial 1 ampl i tude A1 f ven waves (Bj<<B^) are 
linear waves in this sense. The instantaneous Poynting flux for an 
A1 f ven wave i s 


S = c ExB = - -J- (V;xB,)x(B„xB. ) , 
4 TT 4 7T a-'I ✓vO ^ '✓s/0 ^ 

where V, is given by (3). S can be written as 


(B 


li' 


B "B. B, *i*B B,, 
0 111 ~U ■ 0 



4ir 


(5) 


The last term may be written as represents a con* 

vection of the unperturbed field energy by the perturbation velocity 
normal to . For waves which are periodic, the last two terms in (5) 
average to zero, and the average Poynting flux in a wave packet is 


<S> = - 


B, 


Since cose, where = B/ /4 ttp and e is the angle between 

and B, one can write (6). in the form 


4 IT 


* Alfven waves can be classified on the basis of the number of dimen- 
sions in which _B^ fluctuates at an observing point, i.e., there are 
three-dimensional, two-dimensional, and one-dimensional Alfven wave per- 
turbations. Let us consider each type in turn. 

1. Three-dimensional perturbations (General Alfvfen Wave) In addition 
to properties 1.-4. above, a general Alfven .wave has the character! sti cs 

a. l^i = B = 1^0 + JB J == constant. 

b, B, moves in three dimensions. 


.1b1 "CONST 


^.^.^^hcoNST g move on the surface 

/ ^ sphere (Figure la). That this 

n 1 / J Alfven wave can exist has 

/' J/ X y been demonstrated theoretically by 

Whang ( 1973), Goldstein et al . 

"general" transverse (1974), and Barnes (1974). The 

ALFVEN WAVE ALFVEN WAVE „ . . . . 

Poynting flux for this wave is given 
[ by Equation (5). 

X / ' Two-dimensional perturbations 

/ Bo , \ (Transverse Alfven Waves and Rota- ~ 

^y_ j )— tional Discontinuities) If ^]^(t) 

/ fluctuates in a plane, the locus of 

y— y the tip of B, is a circle given by 

SMALL AMPLITUDE ALFVENIC . j- . r- j. [. • T ■ j. L 

ALFVEN WAVE FLUCTUATIONS the intersectioo of this plane with 

Fig. 1. Schematic depiction of the the sphere [B] = constant (Figure 


Thus, for this wave 0, and the 


GENERAL 

alfven wave. 


TRANSVERSE 
ALFVEN WAVE 


SMALL AMPLITUDE 
ALFVEN WAVE 


ALFVENIC 

FLUCTUATIONS 


basic types of Alfven waves (see 
Text). As shown, observed Alfvenic 
fluctuations combine features of 
both general and transverse fluc- 
tuations. 


lb). If the plane does not contain 
the center of the sphere, B^ may be 
defined as the direction normal to 
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the plane. Such Alfven waves are called transverse Alfven waves (TAW's) 
because B, is transverse to B. (B, = 0). If the plane defined by B, 

does contain the center of the sphere 1^| = constant, then B^ moves on 

the arc of a great circle of the sphere, jjB ^ 1 = 1 ^|, and = 0. If 

this case is regarded as the limiting case B^ 0, than ->■ 0 and the 

Poyntin^ flux goes to zero. This means that the "wave" does not propa- 
gate al ong Bq . 

For a TAW there is a minimum variance direction m of £j(t), which is 
parallel to B^ and thus defines B^. Clearly, B^ = ^ • m and, of course, 
B, = constant ^ 0 and B,. = 0. The minimum variance direction is deter- 

mined by finding the direction m which minimizes <[(^ - <£>)■'"] > (see 
Sonnerup, 1971; Burlaga and Turne?, 1976). This is equivalent to diago- 
nalizing the matrix <B.B.> - <B-> <B->, i,j, = 1,2,3. The smallest 

1 J 1 j ^ 2 

eigenvalue, is equal to the variance of ^ along m, i.e., = <Bju > 

A rotational discontinuity is a special kind of transverse Alfven 
wave with the characteristic that B, changes disconti nously in the MHD 

A A ^ 

approximation. The wave moves along B^ = m, which is at an angle e = sifi 
(Bi/Ba) with respect to the ambient field 

3. One-dimensional perturbations (small amplitude waves) If ^ moves on 
an infinitesimal arc of the circle |B| = constant, is parallel to the 
arc and is essentially along one dimension, say p. Clearly this requires 

A 

that Bi << B and it implies that p is normal to <B>. One may choose 

i ^ 

B = <B>, in which case <B.> = 0 (figure Ic). The average Poynting flux 

0 \J n ip A 

is <^> = <B>. This wave has frequently been discussed in the lite- 

rature of small amplitude MHD waves. 

3. OBSERVED CHARACTERISTICS OF ALFVENIC FLUCTUATIONS 

To the Extent that the Alfvenic fluctuations observed in the solar wind 
resemble one or more of the Alfven waves defined above, one can speak 
of the following wave characteristics: 

(a) the correlation between V^ and B^; (b) the polarity; (c) the wave 
speed (d) the locus (dimensionality) of JB^(t); and (e) the wave 

amplitude B^^/B^. Two characteri sti cs of AF's which are not consistent 
W-ith the definition of an Al fvfen wave are: (f) the existence of fluc- 
tuations in , and (g) the presence of tangential discontinuities 
among AF's. It is also commonly assumed that AF's are plane waves, but 
it will be shown in section 4 that this is a questionable assumption. 

We shall now consider in turn (a) through (g) above, 
a. Correlation between and^j^. AF's have been studied in the time 
domai n using both plasma and magnetic field data by Belcher, Davis and 
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Smith (1969) (BDS), Belcher and Davis (1971) (BD), Burlaga and Turner 
(1976) (BT), and Denskat and Burlaga (1977) (DB). Considering 6-hour 
intervals in the period June 14-Nov.21, 1967, BDS and BD found a corre- 
lation coefficient p between and B^^ > 0.8 one third ot the time. 

BD found |p| > 0.6 55X of the time in that period. DB and BT considered 
one-hour intervals in periods during 1967, 1968 and 1971, and they found 
|p| > 0.6 40* of the time. Thus, AF's represent a substantial fraction of 
interplanetary microscale fluctuations. 

AF's have also been identified by Coleman (1966) using cross spectral 
analysis techniques. Considering fluctuations in the range 1 to 50 cycles 
per day for four 10-day periods and two 6-day periods, he found coher- 
ences between and B^ of from 0.20 to 0.70, with an average of 0.46. 

(b) . Pol ari ty The sign of the correlation coefficient or coherence 
between and JBj indicates the direction of propagation of the wave. 

All of the authors quoted in (a) found that essentially all AF's are 
propagating away from the sun. 

(c) . Wave-speed The constant of proportionality between V, and B, is 

K « for an Alfv&n wave. Thus, one can compare the "observed" 

with the "theoretical" value given by = B^/ZAirp by examining the 
distribution of B^K/ ( B^/4irp ) = K/4irp = D. This should be unity for an 

Alfven wave, and indeed BDS, BT, and DB 
found <D> “ 1. However, the distribution 
of D is broad, and values of D>2 are 
sometimes found; the reason for this 
is unknown. 

(d) Locus of B, The distribution of B, 
can be expressed as the distribution of 
the eigenvalue ratios X 2 /A 2 , Xg/Xj, x^/x^ 
which are obtained by diagonalizing the 
matrix <B^. Bj >-<B^ ><Bj > . BDS observed 
that typically one eigenvalue was small 
and the other two larger and of compar- 
able size, indicating that the magnetic 
field fluctuations were primarily in a 
plane, but they did not give the eigen- 
value distributions. BD computed eigen- 
values for the entire Mariner 5 data set. 

Fig. 2. Distributions 9 en- they did not separate AF's from 

values of the variance ellip- 
soids of B for hour intervals other fluctuations. BT and DB obtained 

relative to maximum eigenvalue distributions, shown in Figure 2, of 

X,, showing fluctuation close ^ 

to but not entirely within a eigenvalues of AF's. X- is the variance 
plane. 
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7 ^ 

Figure 2 confirms the result of BDS, 

X3 ^ 0. In particular, BT and DB found that 


0.1. Taken literally, this means that 


of ^ along B^, i.e., = <Bju> 

but it also indicates that 

there are non-zero fluctuations along B^ and therefore that the fluctu- 
ations in B are 3-dimensional, i.e., the waves are general A1 f vfen waves. 
Nevertheless, it is true that the waves are nearl y transverse (Figure 
Ic), and an error analysis has not been made to determine to what extent 
X^ is due to the small instrumental uncertainties. Figure 2 also shows 
that AF's are neither linearly polarized (X2/XJ << 1) nor circularly 
polarized (Bj tracing a circle with X2/XJ * 1). 

(e) . Wave ampl i tude As noted by Hollweg (1978), there seems to be no 
distribution of Bj/B^ for AF'j in the literature. However, this ratio is 
probably close to 0.7 + 0.1 b'ased on studies by Chang and Nishida ( 1973) 
and BD, although those authors considered all types of fluctuations to- 
gether, not just AF's. BT found a most probable valye of between 0.3 and 
0.4 for the relative amplitude of the transverse component of AF's. 

Now we turn to two types of AF observations which are not consistent 
with the properties of Alfvbn waves. 

(f) . Fluctuations in |B| BD observed that small fluctuations in |B| 
could be present among AF's. They found that in general <Bj||><0.1 

and they concluded that Alfven waves are dominant over compressive waves. 

' ' ' ' ' BT showed that in 85% of the AF's they 

examined, 6|^| t 0. Their distribution 
of 6B/B (Figure 3) shows that the average 
is 0.06, which is at least six times 
greater than the experimental uncer- 
tainties. Thus, AF's are not pure Alfv&n 
waves; they must include convected struc- 
tures and/or compressive waves. Consis- 
tent with this are the observations of 
Sari and Behannon (1978) and Sari and 
Valley (1976) that there is generally a 
high coherency between fluctuations in 
|£| and along <JB> and those in at least 
one direction normal to <JB> with phase 
of 0° or 180°. 

( g) . Tangential Discontinuities among 
Alfv&n Fluctuations The conditions for 
AF's are necessary but not sufficient 
conditions for Alfv&n waves. There is no 
guarantee that a discontinuity across 



0.09 01 

8B/B 

Fig. 3. Distribution of rela- 
tive field magnitude fluctua- 
tions for events with 1 p| > 0.6 
and 6B/t < 0.2. It is apparent 
that the fluctuations, though 
smal 1 , are not zero . 
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Which Vj«(Vaq/Bq) Bj is an 
RD, for example. DB showed 
examples of TD's among AF's 
(Figure 4); The TD's were 
identified by determining the 
orientation of the disconti- 
nuity surface from the time 
delay in moving from one 
spacecraft to another, and 
then by showing that there was 
no component of ^ along the 


Fig. 4. A Tu oDserveJ fn magnetic field and ^ that surface. Bur- 

plasma data by both Exp. 33 and Exp. 35.. 

Field is given in terms of solar ecliptic 
components and magnitude as solid lines; 

Solar wind speed, density and temperature 
as dashed lines, Cross-correl ation func- 
tions were computed from magnetic field 
data. 

minimum variance analysis of the current sheets associated with the dis 
continuities. The occurrence of convected structures among AF's might 
simply be the result of Alfven waves moving through an ensemble of sta- 
tionary structures. 


laga et al. (1978) have also 
demonstrated that TD's can 
exist among AF's. They iden- 
tified the TD's by means of a 


4. ARE ALFVENIC FLUCTUATIONS PLANE WAVES? 

It is usually assumed that "Alfven waves" in the solar wind are plane 

A 

waves. This means that B, and V, are functions of x-V„t, where x is the 
direction of propagation. The amplitude, B,, and the propagation vector, 
k, of a plane wave are the same everywhere in the (y,z) plane. Thus a 
spacecraft which intercepts a phase plane at one point should measure 
the same amplitude and k as a second spacecraft which intercepts the 
phase plane at any other point. In other words, two spacecraft separated 
by any distance should measure the same ^(t) profiles (except for a 
phase difference) and the same k. Obviously, truly plane waves do not 
exist. The relevant question is "over what distances are AF's approxi- 
mately planar in the solar wind?". This question cannot be answered with 
observations from just one spacecraft. 

Given observations of Alfvenic fluctuations from 2 spacecraft separa- 
ted by a distance L, one can test for planarity of the waves (1) by 
comparing k^, measured locally at spacecraft 1 with measured locally 
at spacecraft 2, or (2) by comparing k measured locally at one space- 

A 

craft with the phase plane normal, n, calculated on the assumption that 
the wave is a plane wave. If the fluctuations are plane waves, one should 
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find that kj = k 2 and that 2 transverse waves one should 

have that k = m, and n can be calculated by a geometrical argument which 
uses the phase lag of 6B (t) between the two spacecraft and the measure- 
ments of solar wind velocity. 
Actually, with 2 spacecraft 
in the ecliptic plane, only 
the projection of n in the 

A 

ecliptic plane, Op, can be 
determined. This must be 
compared with the correspond- 

A 

ing component of k, i.e.. 


N/M 



< (m. 


n ) 


Fig. 5. Distribution of the angle between 
the eigenvector nip, corres pondi ng to the 
minimum variance direction projected into 
the eclipti.c plane, and the wave vector fi 
for propagating plane wave fluctuations. 


20-70 R^. Their result is shown in Figure 

n was not parallel 

r 


with kp, in testing for 
pi anari ty , 

This latter method was 
used by DB. They tested for 
the planarity of AF's ob- 
served in 1-hour intervals 
by 2 spacecraft separated by 
5. This illustrates that, in 
to kp. Thus the AF's observed 
20-70 


R 


£• 


Since DB 


general , for 17 cases 

were apparently not plane over distances of 
able to test 17 one-hour intervals, it is not known whether their 
are truly representative of all AF's in the solar wind. 


were only 
resul t 


5. PROPAGATION OF ALFVENIC FLUU'UATIONS 


Assuming nearly plane waves and using geometric optics, Barnes (1969) 

A 

and Volk and Alpers (1973) predicted that the wave vector k for MHD 
waves would be refracted as the waves propagate such that k is essen- 
tially radial at 1 AU if it has started out parallel to B„ near the sun. 
This is in the absence of velocity gradients, however, i.e., it is based 
on a simple model of symmetric radial flow. When the velocity gradients 
in streams are also considered, as has been done by Barnes and Hollweg 
(1974), Hollweg (1975), Richter and Olbers (1974) and Richter (1975), 

A 

then, rather than being radially directed at 1 AU, the wave v6;ctor k is 
refracted westward (away from B^) at the front of a stream and, eastward 
(toward B^) at the rear of the stream. 

For plane transverse Alfv&n waves, the wave normal is given by the 
minimum variance direction. Analysis of AF's show that the minimum vari- 
ance direction is usually close to the mean field direction rather than 
in the radial direction (BDS, BD, DB). Solodyna and Belcher (1976), 
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Burlaga and Turner (1976) and Daily (1973) showed that this is true also 
for both compression and rarefaction regions of streams. Thus, the ob- 
served minimum variance directions do not agree with the directions of 
wave normals predicted by the theories. 


FASIA 



PLANE 

WAVES 


SLOW 


AV = 0 
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The theoretical predictions are 
all based on the assumption that the 
waves are plane over distances on 
the order of'^1 AUi as illustrated 
in Figure 6. The tendency for k to 
become radial in the absence of grad- 
ients is the geometrical effect shown 
in Figure 6a. When streams are in- 
cluded, the gradients produce the 
deflections illustrated by the sketch 
in Figure 6b, mainly by convection. 
The various arguments do not apply, 
however, if the waves are not plane 
waves over distances on the order of 
the scale of streams and stream grad- 
ients. In particular, they are inap- 
plicable if the waves are curved 
appreciably on a scale of < 0.01 AU , 
as shown in Figure 6c. For waves that 
are not plane waves, can always 

be close to <B> and still be variable 
over, distances > 0.01 AU because <B> itself is variable over such dis- 
tances. The fact that the autocorrelation length of ^ iss^ 0.01 AU is 
consistent with this picture. 



NON-PLANAR 

WAVES 


SLOW 


(C ) 


Fig. 6. Illustration of (a) plane 
wave propagation without gradients, 
(b) plane wave propagation with 
stream effects included, and (c) 
propagation of waves appreciably 
curved on a scale ^ 0.01 AU. 


6. STABILITY OF ALFVENIC FLUCTUATIONS 

? 

The stability of Alfven waves in the limit g e nkT/(B /8ir)->-0 has been 
investigated theoretically by several authors (Galeev and Oraevs ki i , 

1969; Hollweg, 1971; Cohen and Dewar, 1974; Cohen, 1975; Chin and 
Wentzel , 1972; Lashmore-Da vi es , 1976). Although several types of insta- 
bilities have been identified, the assumption B->0 is not applicable to 
the solar wind, so it is not certain that those instabilities can occur 
i n the sol ar wi nd . 

Recently Goldstein (1978) has investigated the stability of a finite- 
amplitude, circularly-polarized Alfven wave for the case g^l, which is 
the relevant case for the solar wind. He found that the waves can decay 
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for all {5<1 by means of an MHO modul ati onal instability. In particular, 
an A1 fven wave decays into a forward-propagati ng magnetic wave, a for- 
ward-propagating density wave, and a backward-propagating magnetic wave. 
The decay waves are not normal modes of the plasma. The decay occurs as 
a result of the Alfven wave coupling to random density and magnetic fluc- 
tuations in the ambient medium. The growth rate of the driven waves in- 
creases with 6|B|/Band decreases with increasing g. It is surprisingly 
large in general. The decay of AF's may be an important factor in heating 
and accelerating the solar wind. 

The growth of A1 fvSn waves as a result of an instability produced by 
the velocity shear in streams has been investigated by Bavassano et al . 
(1978). They have found that short wavelength waves moving in a direction 
close to that of the mean magnetic field can grow in velocity shear rer- 
gions under the conditions observed near 1 AU, and they observed magnetic 
field fluctuations in such regions which are consistent with those pre- 
dicted by the theory. They also have concluded, however, that most of the 
long wavelength Alfv&n waves could not be produced in this way, and hence 
that they probably originate near the sun. 

7. SUMMARY AND CONCLUSIONS 

Results obtained during the last few years have confirmed that fluc- 
tuations in the bulk velocity of the solar wind are correlated with flue- 
tuations in the magnetic field more than one-third of the time and that 
these AlfvSnic fluctuations have several characteristics in common with 
Alfven waves. However, it has also been shown that AF’s are not "pure" 
Alfven waves, because there are fluctuations in magnetic intensity and 
because tangential discontinuities can be found among AF's. Further- 
more, there is evidence that the plane wave assumption used in most 
theoretical studies of Alfv§n wave propagation is not valid over dis- 
tances > 0.01 AU. 

Several important questions have been raised and should guide further 
studies of AF's. These include the following. What is the cause of the 
variations in magnetic field intensity in AF's? Over what distances are 
AF's approximately plane waves, and how does this change with interplane- 
tary conditions? Is there a turbulent cascade of energy or simple mode- 
coupling near 1 AU? Are AF's general Alfven waves? The answers to these 
questions will require accurate, s i mul ta neous- pi asma and magnetic field 
data from at least two spacecraft with separations in the range 
< 50,000 km to > 500,000 km. Also needed is a better theoretical under- 
standing of the behavior of large amplitude Alfv&n waves, especially 
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with respect to their coupling to other modes, their stability, and 
their propagation. A kinetic theory approach to these problems is par- 
ticularly needed. 
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ABSTRACT A new detailed first principle kinetic theory for solar, v/ind 
electrons which is neither fluid' nor exospheric in nature is 
presented. This theory illustrates the global and local properties of 
the solar wind expansion and coulomb collisions that are reflected in 
the observed electron velocity distribution function — such as its 
bifurcation, skewness, and the "differential" temperatures of the 
thermal and supratherraal subpopulations. This theory is also shown to 
make numerous predictions of correlations amongst electron solar wind 
parameters v/hich have been previously reported. 


The extant models of the solar v/ind represent the extremes of 

approximations that can be made using the Knudsen number (K = 

as an ordering parameter: the fluid approximation (K<<1) enjoys a 

certain popularity because of its tractable mathematics, although the 

transport coefficients (such as the thermal conductivity, viscosity, 

etc.) that are used seem dubious at best (c.f., Kopp and Orrall, 1977, 

for a recent review); the exospheric approach (K “) is hampered by 

the non-uniqueness of the solutions which satisfy the Vlasov equation. 

Direct observations (Ogilvie and Scudder, 1978) indicate that reality 

at 1 AU for electrons is in between these two extremes with K 1. 

(0 

This paper discusses the first quantitative attempt to discuss the 
physics of solar wind electrons in the more realistic, intermediate 
Knudsen number regime. 

The electron Knudsen number at 1 AU is comparable to that eStpected in 
the acceleration regime of the solar wind in the lower corona. If we 
can understand what factor (s) control the observed heat flow, it seems 
likely that this understanding could impact the current reexamination 
of the solar wind energy budget in the acceleration regime. V7e begin 
our discussion with a brief overview of in situ observations of 


electrons at 1 AU which provide yardsticks for assessing theoretical 
attempts to understand solar wind electrons. 

The low energy thermal portion of the electron velocity distribution 
function has been called the "core" or "cold" component and 
approximates, with some license a paraboloid in ln(f) vs velocity; the 
suprathermal component (i.e., everything' which is not the core) is 
often called the "halo" or "hot" component. There is some effective 
Maxwellian with a "differential" temperature which approximates these 
populations. The effective "temperature" of the suprathermal 
component is substantially larger than that of the thermal component; 
its "density" is only a few percent of the ion charge density. The 
signature that this distribution is carrying a heat flux is the 
asymmetry of this distribution function. The electrons that are 
carrying the heat are not those as in small Knudsen number transport 
theory where there is a slight skewness in f(v) in the vicinity of one 
thermal speed; these asymetrie® are observed at higher speeds (usually 
above 2-4 (rms) thermal speeds, and ^observed represented by 

a single Maxwellian. 


A brief list of observed interrelationships of solar wind electrons is 

given below. The halo density fraction of the total di;ti^£iity j,s 

nearly independent of large variations in the ambient d^^hi^ity.* This 

number turns out to be on the order of 6j^2% (where we haVe used 

deviations quoted by Feldman et al., 1975). The ratio of^ the 

— — — — 

suprathermal to thermal "temperatures" is remarkably constant with a 
value something like 6+1. The halo fraction of the total density is 
strongly anti-correlated with the local bulk speed (if it is nearly 
the asymptotic bulk speeds, Scudder and Olbert (1978a)). The 
‘differential temperature of the halo population is anti-correlated 
with the bulk speed; and the heat flux is lower in absolute magnitude 
in high speed wind than in the low speed wind (Feldman et al . , 1976 
and Rosenbauer et al . , 1977). 


A further point of some interest concerns the reproducibly bi-model 
appearance of the electron distribution function. If we plot ln(f) vs 
proper frame kinetic energy for particles moving both along and 
opposed to the local magnetic field direction, there is a break 
energy, E*, between approximately linear segments, which empirically 
scales as E* = (Scudder and Olbert, 1978b) . This energy has an 
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intrinsic relation to the properties of Coulomb collision cross 
section and its velocity dependence. 

The radial variations of the core and the halo components have been 
determined from data returned by the Plasma Science Experiment on 
Mariner 10 (Ogilvie and Scudder, 1978). The radial variation of the 
core electron temperature between 0.4 and 0.85 AU was approximated as 
an inverse power law with exponent 0.3jt0.08 and that of the halo 
component was essentially without radial variation. If we extrapolate 
these radial dependences toward the sun we obtain some indication that 
there is a common zone where these particles may have previously been 

members of the same Maxwellian population above and below the break 

* 

point energy, E . This radial position where the two radial profiles 
converge is in the proximity of 2-15R . 



FIG. la 


FIG. lb 


In Figure la we show 
data acquired by the 
NASA/GSFC VES 
instrument on ISEE-1 
in the format (ln(f) 
vs w„) . The 
characteristic low 
energy thermal 
component and the 
characteristic break 
point energy are 
clearly shown; the 
kinetic information in 
this figure is 
immediate: the system 

is not in local thermodynamic equilibrium (LTE) — if it were in LTE 
with the same internal energy it would look in this format as in 
figure lb. Immediately we can state that the Braginskii (1965) (small 
Knudsen number) formulae for transport coefficients are not relevant 
for this system; Fourier's law (Q = -k VT) does not apply in most of 
the relevant regimes of the solar wind. Fourier's law is a gross 
oversimplification of the physics of heat transport in an 
inhomogeneous plasma and only obtains theoretically under extremely 
simplifying assumptions. We must attempt a kinetic type formulation, 
in the presence of Coulomb collision and Kg'^1; we have to compute the 
electron distribution function both in configuration and velocity 







space. We then have to take the appropriate moment of the distribution 
function to determine the flow of heat. Clearly, if we have a theory 
which can replicate the observation of the distribution function, then 
we may have some hope in developing a quantitative discussion of the 
flow of heat in this system. 

The Coulomb interaction between a "test" electron and the plasma is 
often broken up into regimes for which the test electron has an impact 
parameter, b, less than and greater than the shielding debye length, 

^ ^ , (c.f.j. e.g., Cohen ^ r 1950 or Krall and Trivelpiece 
(1973)). The interactions with b < may be visualized as binary in 
character with shielded ion scattering sites? when b >1^ the 
interaction is a many body interaction of the test particle with the 
correlations of the targets on a scale L > ^ and is a "wave-particle" 
interaction. 

In this paper we will discuss the strict consequences of (binary) 
Coulomb collisions in the cavity occupied by the solar wind. Our 
considerations will give rise to a global sense of causality, v?hereby 
we shall come to understand that the ^ situ electron observations are 
connected in rather dramatic ways to the extremities of the 
heliospheric cavity. This, in turn, implies that to interpret the 
local observation of the electrons in the medium, one must be thinking 
about the system in the large. 

Our approach is to start from the Boltzmann equation with a relaxation 
time approximation for the Coulomb scattering operator. The relaxa- 
tion time approximation is also known as the Krook approximation 
(Gross and Krook (1956)) and is currently being used as a mathematical 
simplification to see qualitative effects? the information obtained is 
already sufficiently instructive to illustrate interesting features of 
a fully ionized plasma which is inhomogeneous. The relaxation time 
scale we have used is the velocity dependent Coulomb time scale for an 
effective 90° deflection; the explicit formulae are discussed in 
detail in Scudder and Olbert (1978a) . The approach is to discuss an 
approach via perturbation where a guess motivated from observation, 
for the equilibrium spatial profile, P(£,v) , is assumed to be a 
convected Maxwellian with a specified radial dependence of density, 
thermal (core) temperature and bulk speed. We then inquire, via the 
relaxation approach, "how do the local Maxwellian populations 


'accommodate' or take notice of their environment • as allowed by the 
smooth forces and Coulomb scattering that the electrons experience, 
enroute to setting up a microscale consistent with the macroscale 
forces and the properties of Coulomb collisions?" The empirically 
motivated choice of a local Maxwellian stems from the evidence at all 
radii for which direct observations have been made that there is an 
approximately Maxwellian core or cold component for the electron 
distribution function. Direct iji situ determinations (Ogilvie and 
Scudder , (1978)) of the mean free path for Coulomb collisions over the 
scale height indicate that there are marginally enough Coulomb 
collisions to start the process of local Mawellization. 


Formal solutions to the Boltzmann partial differential equation are 
rarely of much help in calculating the solution, but they often 
provide valuable qualitative insights into the physical processes at 
hand. The formal steady state solution for '©ur problem at a given 
time, t , at a given location, r , and at a given velocity v, is given 
by 


* * *^ 


* * * 


f<v,r,t )=F(vQ,r^)expC-S^)-t^ exp C-S Cr ^ v; v , t /xlt'] dt' 

t' . 


CD 


where 


SCrivlr*vtt*) = ^ l/xlf’jdt". 


The structure of this solution is particularly interesting since it. 
illustrates the global nature of the eventual solutions we will 
obtain. The first term represents an attenuated, exp(-S^) , signal of 
"source" phase density from the boundary; the additional term 
represents the superposition of "source" phase density which can reach 
the observer from intervening spatial layers properly attenuated for 
the intervening depth, S(r',r ). Let's examine the nature of this 
quantity S in- the formal solution. The quantity S is a counter of the 
effective number of Coulomb collisions experienced by a representative 
element in the phase space going along the characteristic in the 
presence of the smooth forces from (r',t') to (r*,t*) . The quantity 
6xp(“S) which occurs in the formal solution is also a very interesting 
quantity to understand. This quantity is the normalized survival 
probability to transit from (r',f) to (r »t ) without experiencing an 
effective 90° coulomb scattering. The equations of motion are hidden 
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in this formal solution; the attenuation provided by the Coulomb 
effect is explicitly present in the S quantity. The extinction 
provided by the Coulomb scattering is analogous to the similar 
behavior of the optical depth {t{v)) of radiative transfer theory; we 
will refer to this quantity S as the "collisional depth". 

If we consider an observer at lAU it is important to know how large S 
is as a function of energy. This quantity can be shown a posteriori 
to range from approximately 2 to over 30 when one considers electrons 
with local kinetic energy at lAU which range from SOOeV to 5eV. This 
iimnediately implies that in the solar wind we are connected in a. 
dramatic way at higher energies to the plasma microstate in vicinities 
well removed from that of the place' of observation, since the observer 
is not screened off from large phase densities at this boundary at 
this energy. This implies that for any physical quantity which depends 
on the moderate to higher energy electrons, such as the heat flux, one 
must consider the entire system. Because of this global connectivity 
local Knudsen number arguments appear inadequate to vindicate 
Fourier's law in inhomogeneous plasmas. This is not the whole story 
— this is just the extinc'tion portion of the contributions to the 
local phase density; there remains the question of the energetic 
accessibility. 



We have seen that this 

quantity exp(-S) is a 

normalized survival 

probability. Can we understand 

microscopically from the 

properties of Coulomb 

collisions , what allows this 

"transparency" of the corona 

to the observer? In Figure 2 

we illustrate a universal plot 

of the properties of a fully 

ionized plasma for which T « 

e 

(mp/m^)Tp. The ordinate is 
exp(-S) and the abscissa is 
the scaled quantity u = 
w/^T^g, where w is the 
electron's proper frame speed 
and the normalizing factor for 
u is the rms thermal speed of 
the electrons. We have asked 



■FIG. 2 
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the question "what is the survival probability for traversing a 
fiducial local length in the plasma, namely the local free path for 
the rms thermal electron, as a function of test electron speed?" 
Because u is a scaled dimensionless quantity, this is an universal 
plot for the above stated temperature range. We now ask what is a 
reasonable statistical definition that a test electron "survives" 
without having an effective 90® scattering while transiting this 
fiducial length. Since this is a statistical question a customary 

level of confidence is 3 sigma. If we adopt this level of certainty 

* 

we define a characteristic scaled speed u =2.16. If we also recognize 
that this is a statistical relation, then we should expect some 
dispersion in this quantity, say, ranging from 2-4 sigmas. For 
particles above this range of characteristic speeds, we thus argue 
that these particles last Coulomb scattered more than one local rms 
thermal speed electron's free path away from the observer. In 
principle we have established that there is a collisional window above 
this speed range to whatever plasma populations are beyond our 
fiducial length. If we examine breakpoint data from experimental data 
from all the Mariner 10 data and an analysis of published IMP data we 
obtain the most probable points for u* as indicated on the figure. 

These most probable experimental points nearly coincide with the 

★ 

quantity u . The entire range of experimental values of the 
breakpoint energy determined on Mariner 10 is indicated by the 
stippled zone on this figure. This zone is centered in the 2-4 sigma 
zone defined previously. This spread is also consistent with the 
uncertainties in a re-analysis of published IMP data (Scudder and 
Olbert, 1978b). This is a very strong argument in our opinion that 
Coulomb processes may be the dominant factor in determining the local 
distribution function. This also means that Coulomb physics is also 
I intrinsic to the conduction problem since the break point energy is in 
■the proximity of the largest skews of the local distribution function. 

If we are going to integrate the Boltzmann equation we must think 
seriously about the smooth forces' which the electrons experience 
enroute to the observer. This raises the question of what are the 
smooth forces which the electrons experience? Among the most ■ 

important smooth forces is the parallel electric field which the 
charged particles experience because the medium is inhomogeneous. 

This parallel electric field is present in all fluid theories which 
impose the external condition that the medium be pointwise neutral and 
that the electron and ions have the same bulk speeds. This parallel 




electric field can be determined from the generalized Ohm's law (c.f. 
Scuddcr and Olber't, 1978a for details), once the electron pressure and 
density variation have been given (as in our initial ansatz, supra) . 
Our polarization potential barrier is very similar to that published 
by Lemaire and Scherer, 1973). This potential represents a 
significant barrier of the order of 1.5 keV to the progress of 
electrons along the tubes of force from the corona (we have assumed a 
Parker spiral) . This is especially critical since the electrons 
observed at earth presumably originate in a thermal plasma in the 
corona with a kT of the order of lOOeV. 

If we return to the formalism, we can see by integration by parts that 
(1) may be rewritten as a superposition of two terms 

F as - s^) 

fCt} = FCAS = 0) + ^ exp C-&S (Fit CAS)J),)dF C2) 

F CAS = Q) 

CLocal} CGloball 

one that is by definition local, (As=0) , plus a superpositon of 
differential phase density that is accessible to the observer 
multiplied by the probability that it does not get extincted enroute 
to the observer. We see immediately that the departure from the 
local distribution function is this global integral. Since there 
appears, by direct observation, to be a marginally adequate number of 
collisions to Maxwellize the low energy "core" we see from (2) that 
taking F to be the local Maxwellian is a good idea. Among the details 
required to "realize" a solution of the Boltzmann problem, one must 
collect all contributions (in steady state)' which could have ever 
made it to the observer. In this connection it is important to notice 
that in the entire solar system it is very difficult to get electrons 
to cross magnetic tubes of force, either by finite radii effects or 
collisional diffusion. (This means the life and death of the electron 
is along the same field line upon which it was injected in the lower 
corona. This- means that the problem is topologically one 
dimensional). This implies that the stochastic problem of getting all 
the contributions is restricted to the given tube of force which the 
threads the position of the observer. The contributions from a given 
layer at r’ can contribute particles going both fore (away from sun) 

1 and aft ) along the local magnetic field at the 

observer; conversely the observer measuring f^(v ,r ) will in general 
have non-negliglble contributions from radial positions along the tube 



of force that are both inside and outside that of the observer. It is 
quite possible on the basis of kinematic accessibility that the 
electrons which comprise the local aft halo phase density, f , 
originated in the corona, propagated along the local tube of force, 
passed the location of the observer at some earlier time, proceeded 
along the tube of force and have been scattered back toward the site 
where they come to be observed going toward the sun. This same level 
in the corona can also contribute strongly to the forward portion of 
the halo f . Clearly, in steady state there is also the kinematic 
possibility that there are more than one turning point orbits between 
"source" layers and the observer; the problem of getting all the 
contributions comes down to tracking down all these turning point 
orbits. Due to the limitations of space in the present monograph we 
will not pursue this details as they are available in the literature 
(Scudder and Olbert, 1978a). 
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Another example of the 
electron distribution function 
is given in Figure 3 where we 
have plotted the distribution 
function in the proper frame 
adapted from Feldman et al . , 

1975, as discussed in Scudder 
and Olbert (1978a) . The 
principle reason for this 
figure is to show that 
different experimental groups 
obtain similar velocity 
distributions (compare figure 
[la) and establish a basis for 
' comparison with the 
theoretical realizations which 
we will subsequently discuss. 

The thermal ahd suprathermal 
populations are clearly 
present; the break between the 
two parabolas occurs at 
approximately 2 rnis thermal 

speeds; the break in the profiles occurs roughly two and one-half 
orders of magnitude below the maximum value of the entire distribution 
function; the distribution is skewed with a magnitude which should be 
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compared with the theoretical solutions; the ratio of the fore and aft 
phase density f^/f asymptotes to a value of approximately 6. 


Figure 4 shows a plot in the 
same format as that of Figure 
3 of a theoretical solution 
with an asymptotic solar wind 
speed of 400km/s. The 
logarithmically spaced dots 
are the theoretical solution; 
the parabola is the local 
Maxwellian ansatz, F = 

Among the features of this 
profile are the suprathermal 
population, present even 
though no wave particle 
interactions were included in 
this theory, the suprathermal 
population with approximately 
the correct relative 
importance in the overall 
distribution function, the 
skewness of the distribution 
function, which compares 
favorably with that 



FIG. 4 


illustrated in the previous figure of in situ data, the bifurcation of 
the distribution function, and the break- point energy of the solution 
which automatically has the correct position since this is an 
intrinsic property of the Coulomb scattering process. 


From this solution we will establish contact with the experimental 
extrapolation of the radial variation of the core and halo being 
consistent with a prior contact at different times with the same 
Maxwellian population somewhere in the regime of 2^15 R^. Let's ask, 
"where were the collected electrons last members of a Maxwellian 
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'source' population?" We 
consider in Figure 5 these 
radial source regions. The 
vertical axis is the 
contribution of phase density 
which is ultimately collected 
at the indicated observer's 
position and the horizontal 


'SOURCE* REGIONS 


axis is the radial location of 

=t 

the Maxwellian from which the t 


contribution comes. 


The 5eV i 

a: 


seen to have come most 
recently from a Maxwellian 
population in the immediate 


proximity of the observer; 
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f (5eV) comes predominately 
from a Maxwellian population 

FIG 5 

inside the radial position of 
the observer; f~(5eV) comes 
from slightly outside the 

observer. Both these contributions are from within one rms thermal 
free path of the observer. By contrast, consider the radial positions 
of the dominant contribution of the phase density for kinetic energy 
-500eV at lAU. We see that both and f contributions at lAU 
originate in collisional populations very near the base of the corona 


in the region of l-lOB^. We thus see that we recover (by considering 


the Coulomb extinction and the self consistent forces) the 
experimentally inferred regime of common thermal origin, since the 
thermal core has been reprocessed many times in different collisional 
populations with different T's. There is little doubt that at an 
earlier time the thermal electrons were in collisional equilibrium 
with other electrons near the base of the corona. 


This figure also dramatically illustrates the direct contact which the 
lAU observer has (at suprathermal kinetic energies) to the boundaries 
of the system. In these senses the ^ situ electron distribution 
function has both local and global information in it: the thermal 
regime has the local gradient physics; the population of the f"*” in the 
suprathermal regime proclaims the existence of the hot, very dense 




corona many mean free paths away- and the skewness of the 
super athermals contains the information of the transit of these 
particles beyond the observer in the distant heliosphere where they 
are backscattered to be collected going toward the sun. 

The typical radial position where these particles are turned to move 
again toward the sun is 5-lOAU, depending on the energy. We think the 
key point that has been missed in the exospheric calculations is that 
the wrapped nature of the interplanetary field in the ecliptic 
introduces a longer arc length path to exit the heliospheric cavity 
than estimates based on the characteristic radial dimension of the 
cavity. This in turn, implies that the Coulomb mean free path and the 
free path of any energy electron divided by the scale height along the 
magnetic field decreases with increasing distance; this in turn 
implies that it is inconsistent to neglect the backscatter of 
energetically unbound electron back into the heliospheric cavity. 

(Cf . Scudder and Olbert 1978b and Sittler et al. (1978) (this 
monograph) for further discussion on this point.) 


To establish some contact with 
the parameter variation with 
the asymptotic solar wind 
speed which we noted in the 
introduction, we consider in 
Figure 6 a theoretical 
solution for an asymptotic 
solar wind speed of 800km/s. 

It can be shown in a rather 
general way that faster wind 
implies a higher electrostatic 
potential energy barrier for 
the coronal electrons to get 
to the lAU observer. It is 
thus not surprising to see 
that the suprathermal 
populations are depleted 
relative to the thermal phase 
density-this was one of the 
strong anti-correlations which 
has been noted in the data. 

If one compares this solution 
with that of Figure 4 we also 



FIG. 6 


that the effective suprathermal 
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"differential" temperature is lower in high speed wind than high. 

This is also a reported correlation in the solar wind electrons. The 
lower mean energy or differential "temperature" for the suprathermals 
in the high speed wind arises principally as a result of the dominant 
radial position of ballistic connect'ion shifting into the upper corona 
(10-20R ) as a result of the enhanced barrier between the lower corona 
and the observer. Since the heat flux is an average over three 
properties of the distribution functions Q = ®'^skew^skew ’ 

(where n is the number of particles carrying the heat, T is their 
effective mean energy, skew is the asymmetry of the heat carrying 
regime of f (v) , and a is a dimensionless constant of order 1) , and all 
of these quantities are inversely correlated with the asymptotic bulk 
speed we immediately have that the local heat flux should be 
anti-correlated with the asymptotic bulk speed. This is a global 
relation and has none of the proximate character implied in the 
Fourier's heat law that is so commonly used. This latter . 
anti-correlation has been reported at lAU by Feldman ^ (1976) in 

another context and also reported by Rosenbauer £t (1977) inside 
lAU. 


For lack of space in this monograph the details of the empirical 
corroborations of the expectations of this theory may be found in a 
forthcoming series of papers in the JGR: (Scudder and Olbert 1978a, 
1978b). 

CONCLUSIONS ; 

Fourier's law is certainly not a correct ■ description for the transport 
of heat in an inhomogeneous plasma. The thermally driven stellar 
(solar) wind problem should be reposed. The principle impact will be 
for the acceleration region where perhaps inadequate estimates have 
been previously made for the rate at which internal energy can be 
supplied to the expanding flow. Coulomb collisions play an essential 
role in the transport physics of solar wind electrons. 


The observation of a suprathermal population in an inhomogeneous 
plasma cannot be taken as a priori evidence of wave-particle 
interactions; by the nature of our construction we have excluded wave 
particle interactions, yet the combination of inhomogeneity and the 
Coulomb "window" above E*=7kTg naturally gives rise to this leakage of 
non-local collisional populations in superposition with local 




collisional populations. This is at odds with the Parker and Tidman 
(1958) suggestion tihat all suprathermal populations are generated as a 
result of Fermi or betatron action. 


A theory necessary to show the global signatures in the _in situ 
electron data is in its formative stages; these techniques have 
already suggested new ways to look at the data to extract information 
about the inaccessible regions of the solar wind expansion. 


The theory that has been developed makes definite predictions, which 
have been corroborated. As an example the theory predicts (Scudder 
and Olbert, 1978b) the asymptotic radial variation for the electrons 
thermal temperature to be inverse 1/3 power as is supported by the 
Voyager 2 electron observations reported in this monograph by Sittler 
et , (1978) . 
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ABSTRACT Preliminary estimates of the electron core temperature as a 
function of radial distance from the sun between 1,36 and 2.25 AU are 
presented. The electron core temperature is found to gradually 
decrease with a power law index equal to -0.34 + 0.16. 

INSTRUMENT 


The Voyager spacecraft plasma instruments measure solar wind electrons 
in addition to the ions as have been reported by Belcher (1978) 

in this monograph. The electron sensor (see Figure 1 in Belcher £t 
al. , 1978) is a cylindrical, potential modulated Faraday cup, with 
aperture normal perpendicular to the radial and 23° above the ecliptic 
plane. The field of view of the sensor is conical, with a half angle 
of 45°. The angle between the average Parker magnetic field direction 
and this normal is a decreasing function of the radial distance, with 
a maximum value of 45° near earth and a minimum value of 23° at 2.25 
AU as the interplanetary m.agnetic field becomes increasingly 
transverse to the radial. The electron measurements are derived from 
fluxes which are locally directed av/ay from the sun along magnetic 
tubes of force. 


The electron sensors have tv/o energy modes El: 10 eV - 140 eV and E2: 

140 eV - 5950 eV, each mode possessing 16 contiguously spaced energy 

channels with nearly equal widths of 320 km/s in speed. We will 

concentrate on the low energy regime in which the thermal component of 

the electrons are usually found. These are the first electron 

measurements which are differential, yet contiguous across the thermal 

plasma regime. This narrow spacing allows a determination of electron 

4 o 

temperatures as low as 3 x 10 K, where the electron thermal speed is 
1000 km/sec. 
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Two samples of electron 
spectra from Voyager 2 
data, acquired at 1.47 
and 2.38 AU, 
respectively, are shown 
in Figure 1. These 
spectra were selected as 
they represent similar 
flow states of the solar 
wind at different radial 
distances. The actual 
macroscopic parameters 
are given beneath each 
figure. The ordinate of 
each spectrum is the 
reduced distribution 
function, F^(v*n), given 
by 
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where n is the unit normal of the sensor, v^ is the mean tuned speed, 
and A is a proportionality constant. The horizontal axis of each 
figure is the apparent speed of. the detected electrons at the 
spacecraft. This speed is related to the ambient electron's speed 1 v^| 
before encountering the spacecraft sheath by the relation 


where 4> is the spacecraft's electrostatic potential with respect to 
the plasma. The solid line histogram form of the plotted data 
explicitly shows that the measured currents, represent integrals 

over a speed range, 


Since the current, I, collected by the detector represents a folding 
of instrument characteristics with the ambient properties, the reduced 
velocity distribution, is Tiot strictly proportional to the 

ambient phase density, It can be shown that F^(Vj^) is very 

similar to the actual velocity distribution. This is best illustrated 

by considering the crosses overlayed on this figure. The crosses 

represent the best fit convected Maxwellian, in the presence of a 

spacecraft potential, convolved with the instrument's response function 

for the thermal electrons (core) . A Maxwellian in this format would .o 

ORIGINAL PAGE IS 

OP POOR QUALITY 





41 


be a parabola. The convolution of the Maxwellian is almost 

indistinguishable from, a parabola. We are thus justified in 

discussing the shape of as being very nearly that of so 

long as f (v) does not have structures with half angles of less than 
_ ® 

^,40°. In both spectra we see in a characteristic low energy 

regime, where the convolved thermal Maxwellian (+'s) agree very well 
with the observations; these observations also show a no:n~Maxwellian 
suprathermal tail. For these two examples the departures from the 
thermal Maxwellian occurs at a speed for which the electrons have 
characteristic kinetic energies of 6 - 7kT^ where is the thermal 
temperature and k is Boltzmann's constant. The suprathermal electrons 
are almost always detectable above the one telemetry count threshold 
which itself generally represents an ambient signal larger than 
detector noise levels. In the more distant spectra the detector gain 
has been increased by a factor of 10, yielding a corresponding 
increase in sensitivity per telemefct^y count. These two representative 
examples clearly demonstrate the instruments capability to meaning- 
fully measure and allow analysis of electrons many astronomical units 
from the sun; in this connection, as of this writing, Voyager 1 is 
-acquiring electron data well above one telemetry count and instrument 
noise beyond 4 AU. 

ANALYSIS 

We have determined ambient thermal electron parameters under the 
assumptions that: (1) the thermal electrons are well represented as 

an (isotropic) convected Maxwellian (2) and that the principal effect 
of the spacecraft potential on the sampled distribution function is 
that it is shifted in energy by the relation of equation (2) (cf., 

Sittler, 1978 for discussion of this point). The first assumption 
will be supported a posteriori by determining the local Coulomb mean 
free path over the scale height (Knudsen number) and showing that 
there are an Increasingly adequate number of pitch angle scatterings 
to isotropize the thermal electrons. Because the core electrons 
dominate the electron density they cannot have a common bulk speed 
very different than that of the ions. This has been verified at 1 AU 
by Montgomery et al. (1968) and Feldman et £l. (1975) . The second 
assumption follows from an extension of Mott-Sraith and Langmuir's 
(1927) results in the presence of highly subsonic flows (qf., Sittler, 
1978) . Since some attempts have been made to make the Voyager 
spacecrafts near iso-potentials to prevent arcing at Jupiter 
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encounters, we have an a posteriori check on the reasonable magnitude 
of the derived spacecraft potential. Once we have made these 
assumptions an analytical expression for can be written 

which is a non-linear functional of the core electron temperature, , 
and the spacecraft potential, The electron density is not a free 

parameter in our analysis, since we fixed the normalization of the 
Maxwellian by requiring it to equal the positive ionic charge density 
determined from the forward cluster' (n = n + 2 n ) . We acquired a 

® P -1 “ 

zeroth order estimate of T ' as (d £nF /dE) /k where the 

c e» ' ' 

logarithmic gradient is determined using the second and third speed 

channels. Note, the first energy channel has been temporarily omitted 

from our fits for now, since it has not yet been calibrated properly. 

We then use the first five measured channels (2-6) to determine, via 

gradient search, the optimal and . In general these five 

point fits were poor since the higher channels contained an admixture 

of suprathermal populations. We kept deleting the highest speed 

2 

channel and redoing the f.i.t until it was either x acceptable or we 

were down to three points (one degree of freedom). In some instances 

4 o 

when T ^ 4 x 10 K, the three points fit were unacceptable and were 

Omitted from subsequently averaged data. The reason for this 

difficulty which confined our analysis inside 2.25 AU is not due to 

the instrument's inability to resolve such a steep spectrum, but is 

due to the core electrons beginning to move inside the lowest energy 

channel as they cool. This effect is consistent with the scaling law 

, ★ 

for the breakpoint energy E = 7kT proposed by Scudder and Olbert 

c 

(1978). Because of this, one may have only one or two points 
characteristic of core electrons, while the remaining are 
characteristic of suprathermal electrons. This selection effect, 
which in our opinion does not seriously affect the current radial 
domain, would have caused our daily averages to be biased toward 
higher temperatures, especially at larger radial distances where the 
electrons are getting colder, when the first energy channel does 
become available our analysis will be carried out to larger radial 
distances. 
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An overview of the radial 
variation of thermal 
solar wind electron 
properties is presented 
in Figure 2. We have 
plotted daily averages of 
the Coulomb-Knudsen 
number CAjjj^p/H) , the 
electron "core" 



FIG. 2 


temperature, T , the 

proton number density, 

n , and the solar wind 
P 

(proton) speed. In 
support of our first 
assumption above we see 
that the Knudsen number 
is less than one for all 
stream conditions and in 
the radial interval has 

an average of 0.3. The scale height that has been used for the 
plotted Knudsen number is H = (d i,n n^/dr) ^ is actually a lower 
bound, since the more appropriate scale height should be H|j = (d 
Hg/d ^ where the derivative is taken along the local magnetic field 
direction. These two facts imply that on average in this range the 
thermal electrons experience in excess of 4 collisions per scale 
length, Hj| , which is sufficient on theoretical grounds (Spitzer, 

1953) to "Maxwellize" the core electrons. Although not shown for 
lack of space, the inferred, median spacecraft-plasma potential is 5.7 
eV, consistent with theoretical estimates of this quantity for a 
conducting spacecraft surface in sunlight. 


The core temperature shows the usual compressive heating in stream 

interaction regions and the tendency to be locally cooler in high 

speed than low, that would* be expected for its behavior like a 

classical gas (cf.> Scudder and Olbert, 1978a, b) . The electron core 

5 o 4 o 

temperatures vary between 1.5 x 10 K and 4.5 x 10 K with an average 

4 o 

value of 7.4 x 10 K, that is significantly lower than the mean value 
of 1.25 + 0.29 X 10^°K reported by Feldman £t al. (1975) based on 1 AU 
observations. 
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The daily average temperature data indicate as a whole a gradual 

decrease with radius, best fit (with formal errors) by an inverse . 

power law with index a ~ 0.34 + 0.16 and 1 AU intercept values of 

(1 AU) - 10^*^^ In an effort to try to reduce stream 

interaction effects we also computed sliding linear least square fits 

to 3 data points which were spaced 26 days apart, none of which being 

within- a leading edge of a stream. We thus obtained a distribution of 

slopes and intercepts giving a weighted mean inverse power' law index 

4o 

of 0.47 and temperature at 1 AO equal to 9.S :< 10 K. These values 
happen to be in good agreement with our average fit to all the data. 

A 

We also fit only the first two solar rotations of data to the power 
law relation to minimize the selection effect discussed previously. 

We obtained a power law index consistent with that reported for the 
entire three rotation fit. 

A full analysis of this data set has not been completed. In order to 
present a timely preliminary survey for the Solar Wind 4 Conference 
all this data have been processed in the last two months. In time, a 
more thorough evaluation will be possible. 


DISCUSSION 

The preliminary radial variation of the thermal electrons between 1.36 

and 2.25 AU of r , a = 0.34 + 0.16 is consistent v/ith the recent 

determination of this quantity by Ogilvie and Scudder (1978) as 

determined from Mariner 10 data acquired between 0.4 and 0.85 AU, and 

as determined by Gringauz and Verigin (1975) between earth and orbit 

of Mars. This asymptotic variation in the core temperature is 

centered in the range (0.21 - 0.44) suggested by Scudder and Olbert 

(1978b) in order that the Coulomb collisional window at E* = 7kT be 

c 

the local polarization potential energy to infinity. This radial 
variation appears to imply that the core is becoming more and more 
collisionless since the gas Knudsen number scales as 

- - S' 1^-2“ 

r n 


V- _ 4T « 1 t2 

S H 


for inverse square density variation and inverse power law index a for 
the temperature variation. If a < 0,5 then K_ grows with radius. The 
. best fit (x= 0.34 seems to imply by this type of argument that the 
core will become less and less collisional with increasing radius. 
However, the magnetic Knudsen number for any latitude not over the 
solar spin axis goes as 

K„ - ^ 

«l! 


which decreases for any inverse temperature variation. We thus 
conclude that our data are consistent with an increasingly collisional 
character of the core electrons as r increases. 
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OBSERVATIONS OF THE *STRAKL» DY THE SOLAR WIMP ELECTRON 
SPECTROMETER ON MARINER 10 


K. W. Ogilvie 
J. D. Scuddec 

NASA/Gcddard Space Flight Center 
Laboratory for Extraterrestrial Physics 
Greenbelt, MD 20771 

ABSTRACT A search of Mariner 10 electron data has been made for 
examples of"unusual distribution functions, and three examples found 
which are-identified wi:th-the-”Stc-ahl" previously described by 
Rosenbauer-et al. (1976). The electron sprectrometer on Mariner 10 
observed only. electrons in motion .towards* the sun so the present 
results may be due to a "kinked" magnetic field geometry. The 
characteristic large increases in electron flux density at energies 
above about 50 eV appeared when the magnetic field direction was 
within the solid angle of the instrument and disappeared when this 
condition was -no longer fulfilled-. At other times and distances, the 
magnetic field direction, was often.within the solid angle of the 
instrument at times when no "Strahl" was observed. All three 
observations were made early in the encounter of the spacecraft with a 
high-speed stream, at heliocentric distances of between 0.65 and 0.45 
AU. The angular response of the instrument is 43.5° in the ecliptic 
plane and +13.5° perpendicular to that plane, and contour plots 
indicate the "Strahl" to be a narrow feature with, an energy dependent 
PWHM ranging between 13 degrees at 400 eV and r- 40° at 50 eV, 
superimposed upon the normally guasi-isotropic halo part of the 
electron distribution function. 

INTRODUCTION 


This paper will describe and interpret observations of a distortion of 
the interplanetary electron distribution function along the direction 
of the magnetic field. This distortion was first described by 
Rosenbauer et al. (1976) , when discussing observations by the plasma 
instruments on the HELIOS spacecraft, and called by him the 'strahl^. 
The present observations were made by the plasma electron spectrometer 
on Mariner 10 (Ogilvie et , 1977) at radial distances between 0.45 
and 0,65 AU. This instrument scanned through an arc of 120°, from 
^ = 140° to = 250°, Figure 1, about the anti-solar direction, 
and thus could only observe electrons in motion towards the sun. A 
search of the. 90 days of data tiaken during th-a cruise from 1 AU to 
Mercury turned up at least three examples of the phenomenon. The 
relationship of this distortion to the magnetic field is illustrated 
by the fact, that the strahl was only detected by Mariner 10 when the 
magnetic field vector passed v;ithin the sensitive cone of the 
instrument, bringing with it the axis of the strahl. The effect was 


not seen, however, every time the magnetic field direction was within, 
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this cone. The examples of the strahl found in the Mariner 10 data 
supported the contention of Rosenbauer e_t (1976) , that the strahl 

is usually found early in the observation of a high speed stream. 

j 

OBSERVATIONS 

The electron plasma instrument on Mariner 10 has been fully described 
by Ogilvie et al . (1977) ; the most important characteristics for the 
present purpose are the energy range (13.6 to 700 eV) and the angular 
resolution (+3,5 degrees in a plane perpendicular to the scan axis and 
+13.5 degrees in a plane containing the scan axis). The scanning 
motion was at a rate- o£ one degree per second about an axis. Figure 1, 
perpendicular to the ecliptic plane. A .15 channel differential 
electron energy spectrum was taken once every six seconds. The 
Mariner 10 data was searched for spectra with abnormally high fluxes 
at high energies, and the events described here are part of the 
results of that 
search. Figure 2 
shows the most 
extensive of the 
five observations, 
and the one from 
Which the most 
information can be 
extracted. The top 
5 panels are 
logarithmic plots of 
electron counting 
rates during a 3 1/2 
hour period, 
measured when 
Mariner 10 was 
situated at a radial 
distance of 0.51 AU. 

At this time the 


TO NORTH 



FIG, 1. The geometry of the detector . 

Electrons with velocity vectors falling inside 
the angular space lei < 13.5°, .0 Q to 70 ' 
degrees and 0 = 32u to zero degrees can enter 
the detector at some time during its angular 
scan. The scan axis points within a degree of 
Korth, 
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. 2 

. . Observation of 

strahl on' March 16, 1975- 

The double 


peaks at the higher electron energies are caused by the scanning 
motion.” The angle 6 is the magnetic field latitude angle.- The 
angle 6 is the angle between the detector and the magnetic field 

spacecraft, had just encountered a compression region,. and was situated 
in a positive magnetic sector. Until about 2130 UT the low energy 
channels (represented by the 40 eV and 71 eV channels in Figure 2) 

..... ■ j * 

showed a small flux variation with scan angle but at higher energies 
(220 and 389 eV, for example) no such variation was present and the 
fluxes were low. Between 2115 and 2140 UT, the interplanetary, 
magnetic field direction rotated slowly from a direction 0 = 236®, 6 = 
30°, to a direction 0 335°, 0 = 0°. Electrons can enter the instrument 
at some time- during its scan if their velocity vectors fall in the 
angular regions 0 = 0 to 70 degrees and 0 = 320 degrees to 0, 0 ^ 6s 

13.5 degrees. Figure 1, so that electrons aligned with the magnetic 
field could then be observed, and greatly increased highly anisotropic 
fluxes appeared in the higher energy channels; there was little change 
at lower energies. The characteristic double-peaked appearance at 
these energies is introduced by the motion of the detector, which 
traversed the velocity space structure tv;ice in each angular scan. In 
response to a slow change in the magnetic field direction, to 0 = 

308°, 6 = -26°, the flux of field aligned electrons again disappeared, 
We can see from Figure 1 that the angular width of the 
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strahl decreases as energy 
increases. Figure 3 shows a 
normarized logarithmic contour 
plot (of Log 
distribution function 
projected onto the ecliptic 
plane for- a time period of 
three minutes at 2130 UT 
during this event. The 
quantity fj^^ is the value of 
the distribution function 
derived from the lowest energy 
channel (13 ..6 eV) . The vector' 
labeled B represents the 
direction of the magnetic 
field, projected on the scan 
plane, and the shaded sector 
shows the angular width (7®) 

Qf the instrument acceptance 
la this plane. One can see 
from Figure 3 that the 
apparent PWHM of the strahl, 
at an energy of 350 eV, is of 
order thirteen degrees, about 
twice the angular resolution 
of the instrument. Correcting 
for the non differential 
instrument response yields a 
value of 12.9° for the FWHM, 
Thus the strahl appears as a 
distortion of the distribution 
function v;hich can be measured 
at all energies but is most 
marked at several hundred 
electron volts, and therefore 
characteristic of the halo. 

It is very strongly correlated 
with the magnetic field 
direction. In Figure 4 we 
plot the distribution function 
along the axis of the feature 
and perpendicular to that 



FI(j, 3, Logsirthralc contour plot of the 
electron distribution, function .for the ' 
event in Figure 2 at 2133 UT, 



PXG, -4, Th,e distribution function 
along th.e, axis of the strahl and 
approximately perpendicular to It, 


It 


9 







(dxis. We see that the 
differential temperature 
characterizing the core 
electrons up to 50 eV shows 
little change, while that of 
the higher energy electrons is 
almost twice as great along 
the axis as perpendicular to 
it. The distribution function 
also appears to be a better 
fit to a Maxwellian along the 
axis than perpendicular to it. 
-la Figure 5 we see another 
example, an observation made 
on February 22, 1974, at a 
heliocentric. distance of 0.635 
AtJ in front of a high speed 
stream. There the magnetic 
field rotates from 0-90 
degrees, 0 “*0 degrees to 0 = 
135 degrees 0 = 0 degrees and 
back again. This is about 
five degrees less than 
required for the magnetic 
field direction to enter the 
instrument, but the 
anisotropic counting rate at 
high energies appears due to 
the finite velocity space 
width of the structure. 
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Table I shows details of the five strahl detections; the search of the 
Mariner 10 observations during which these were found was conducted 
for unusually high counting rates at high energies. We cannot 
therefore state that no other strahl observations exist undetected in 
the data, but' our purpose is not to survey general conditions of 
occurrence of the strahl, for which a sunward-facing instrument is 
obviously more suitable, but to use the angular resolution of the 
present instrument to study a few cases, primarily the March 16, 1974 
observation shown in Figure 2. 



TABLE I 
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Date 

Radial 

Position 

Relation to 
Stream 

Remarks 

1. 

Feb. 22, 1974 

0.635 

Immediately before 

Figure 5 

2. 

March 3-4, 1974 

0.59 

field maximum 
Questionable 


3. 

March 16, 1977 

0.51 

Immediately before 

Figure 2 

4 . 

March 22, 197‘4 

0.45 

field maximum 
Immediately before 


5. 

March 24, 1977 

0.48 

field maximum 
High bulk speed 


CONCLUSIONS- 





The observations show that electron distribution functions having 

distortions along the magnetic field direction can be observed by an 

anti-sunward facing detector. The angular half-width of the strahl 

decreases with increasing energy, with a FWHM value of 12. 9 degrees at 

ah energy of 350 eV, and the electrons measured along the magnetic 

5 o 

field have a temperature of order 8 x 10 K. Such distribution 
functions are observed at positions close to the density increase at 
the leading edge- of high-speed streams, in general agreement with 
Rosenbauer (1976) , In Figure 6 we see two hypothetical geometrical 
arrangements leading to the observation of a strahl in the anti-solar 
direction. Iri each case the spacecraft is situated at the point X, 
just in front of the compression region of a high speed stream 
advancing towards the left. In the first case, a, the strahl is 
formed by electrons moving away from the sun in a kinked field 
configuration. The variation of the magnitude of the magnetic field 
along a field line passing through the point of observation is shown 
schematically in the B-L graph; the locally collisionless halo 
electron component forms the strahl in the region of decreasing 
magnetic field. In case b, there is no kinked field configuration, 
but electrons scattered back towards the sun at great distances are 
focussed into a strahl by moving from the higher magnetic field in the 
compression region into the region of lower magnetic field before it. 
This is again illustrated in the B-L graph inset into Figure 6b. A 
decision between these two possibilities must await the interpretation 
of more examples. 


FIG . 6 , Tvco possible magnetic field conf igU’«^ 
rations wbicii could be associated with th.e’ 
observation of a strabl In the cases described in 
this paper. The inset graphs indicate the vari- 
ation of the magnetic field B along the arc -line 
L, going from A to B. 
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An Evaluation of Corotating Solar V7in<3 Stream Models 

V. Pizzo 

NAS/NRC Resident Research Associate 
Laboratory for Extraterrestrial Physics 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 USA 

ABSTRACT 

Vie conipare the predictions of six nonlinear, single-fluid, corotating 
stream models ranging in sophistication from a one-dimensional 
hydrodynamic description (1-D HD) to the full three-dimensional 
magnetohydrodynamic (3-D MHD) formulation, for initial conditions at 
0.3 AU intended to mimic the Helios close-pass data. It is shown that 
inclusion of both the multi-dimensional flow and the magnetic field 
are crucial to accurate description of the resulting stream evolution 
inside 1.0 AU. For most dynamical studies, the simpler 2-D MHD model 
will suffice. The full 3-D version is required only to investigate 
those properties strongly affected by the north-south flow. 

Two main elements comprise any numerical simulation: first, the 
physical processes implicit in the equations used to describe the 
system; and second, the input conditions. Depending upon the input 
conditions, the relative importance of the physical mechanisms in the 
model can change dramatically. Therefore, before we can safely apply 
a model to the analysis of complex real data, it is necessary to 
engage in a systematic investigation into the properties of the model, 
using carefully chosen hypothetical examples that are both simple 
enough to facilitate analysis, yet realistic enough to touch upon the 
relevant physics of the system under study. 

It is our purpose here to examine a family of basic single-fluid 
models of corotating structure in the solar wind — ranging in 
sophistication from a 1-D HD formulation to a 3~D MHD description — to 
see how the dimensionality of the formulation and the treatment of the 
magnetic field relate to the predictions forthcoming from these 
models. 

Steady, corotating solar wind flow in the supersonic 
interplanetary regime may be approximately described by the MHD 
equations, presented here in vector form in the inertial frame: 

^ r 

4 U*V) = 0 

PAQE BLANK NOT FILMED 
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~ =• Vx(ux3) 


The independent variables are the usual spherical coordinates (r, 

Q f while the dependent variables are the velocity, u, the density, 
p, the pressure, P, and the magnetic field, B. The polytrope index, y , 
is set at 5/3 in the examples presented here, G is the solar 
gravitational constant, and n is the solar equatorial rotation rate. 
With the assumption of field-aligned flow in the rotating framer the 
last (induction) equation may be eliminated entirely (except in the 
Ir-D MHD formulation). Thus, in the full 3-D MHD case, we are faced 
with a coupled system of six hyperbolic, nonlinear, PDE's in three 
special variables. Details of the numerical techniques required to 
solve such systems may be found in Pizzo (1978a, b) , 

The essential physics embodied in these expressions is best viewed 
in the second (momentum) equation. The two terms on the left contain 
the basic kinematic steepening effect known as the stream interaction 
mechanism. (The term explicitly accounts for the exact 

relation between time and azimuth for corotating structure in the 
inertial frame.) The first two terms on the right describe the 
dynamical reaction of the plasma to these rotational-inertial forces. 
Thus, there exists a competition between the kinematic steepening 
mechanism on the one hand, and the dynamical broadening processes 
attributable to gas and magnetic pressure forces (and, to a lesser 
extent, magnetic tension) on the other. 

Within the context of MHD models, the two main factors regulating 
the course of the dynamical evolution are the intensity of the field 
and the dimensionality of the formulation (the number of degrees of 
freedom allov;ed to the fluid motions) . On this basis, we can 
establish a classification scheme encompassing six models, most of 
which, in some variation, have appeared before (1-D HD: Hundhausen, 
1973; 1-D MHD: Steinolfson et al . , 1975; 2-D HD: Goldstein, 1971; 2-D 
MHD: Nakagawa and Wellck, 1973, Goldstein and Jokipii, 1977, and Han, 
1977; 3-D HD; Pizzo, 1978a; 3-D MHD: Riesebieter, 1977, and Pizzo, 

1978 b) . 
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Our procedure will be to compare solutions obtained with each of 
the six models for identical input conditions, beginning at 0.3 AU and - 
terminating at 1.0 AU. Thus the results of this paper apply only to 
stream dynamics inside earth orbit; a follow-on will deal with ' 
corotating structure in the distant solar wind, v/here the streams 

are highly evolved and shocks dominate the flow. 

In. Fig. Ir va show the 
input variations at 0.3 
AO in radial velocity, u^, 
number density,, fi, and 
temperature, T as a 
function of solar 
longitude. The nature of 
this hypothetical stream 
distribution — the narrow 
transition between the 
flow states and the marked 
intrinsic temperature- 
density-correlations — 
were chosen to mimic the 
Helios close-pass data 
(Rosenbauer et al . , 1977). 

The other parameters, as 
appropriate in each model, 
were treated as follows; 
the gas pressure is held 
constant across the stream 
at the inner boundary, as 
is the field intensity (jBl 
= 45y) ; the azimuthal 

velocity, u^, is chosen such that the flow tubes in the rotating frame 
are all initially parallel; the meridional velocity, u is set equal 
to zero; and the geometry of the 3-D distribution is a simple figure of 
revolution centered at the equator (i.e., circular contours on the 
source surface r^ = 0.3 AU) . For numerical convenience, we postulate 
periodic boundary conditions in the azimuthal direction. None of these 
latter parameter specifications seriously affect our results. 
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FIG. 1. SOURCE CONDITIONS AT 0.3 AU 
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Given these initial conditions at 0.3 AUr we use a 
finite-differencing form of the MHD equations to simulate the stream 
(‘volution to 1 AU. Figure 
2 depicts solutions in the 
equatorial plane at 1 AU 
for the 1-D, 2-D, and 3-D 
MHD models. Just two 
parameters, the radial 
velocity and the number 
density, suffice to 
contrast the behavior of 
the models. Ke 
immediately see that the 
3-D and 2-D solutions are 
quite similar, whereas the 
1-D prediction is 
radically different. What 
has happened is that the 
nonrad ial flow in the 
multi-dimensional cases, 
though small (Ug, u^<<u^) , 
acts to smooth out the 
compression significantly. 

As the fluid is radially 
squeezed in the stream 
interaction, a vigorous 
shear flow develops at the 
interface between the fast 
laterally, thereby relieving the compressional stresses. This 
nonradial flow is not included in the 1-D description, and hence the 
evolutionary state of the stream is grossly overestimated. (In this 
case, the 1-D model predicts shock formation at about 0.4 AU, contrary 
to the observations, while the multi-D models shock near 1.2 AU.) 

Figure 2 also demonstrates that the burden of the nonradial flow 

effects are contained within the 2-D formulation; little is gained by 

expanding to 3-D. This behavior arises because the azimuthal flow is 

more directly driven by the rotational interaction than the meridional 

flow, which is of a decidedly subsidiary nature. The north-south 

transport of material spurred by the meridional motions has little 

effect upon the gross stream dynamics at any given latitude, except at 

the very fringes of the stream, where the amplitudes are low. ORIGINAL PAGE IS 

Of- POOR QUALITY 
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Fig. 2 is a direct measure of dimensionality effects in stream 
dynamics. However, an equally important question concerns the role of 
■the smooth interplanetary magnetic field (IMF) . Because the field is 
tied tp the fluid motions and because the thermal and magnetic energy 
densities are comparable in the solar wind, we might expect magnetic 
effects to be significant in the interaction dynamics. 

in Fig. 3, we isolate 
the magnetic influence by 
considering 2-D MHD {solid 
line, same as Fig. 2) and 
2-D HD (dashed) solutions 
at 1 AU executed for the 
input conditions of Fig. 

1. We see that, like the 
1-D MHD model, the 2-D HD 
description overestimates 
the kinematic steepening, 
with a forward-reverse 
shock pair, forming near 
0.8 AU. (The compression 
region in the 2-D model is 
narrower than in the 1-D 
MHD at 1.0 AU since that 
of the latter shocked 
close to the sun and 
expanded thereafter.) The 
presence of the IMF 
approximately doubles the 
local characteristic speed 
at which disturbances 
propagate, thus spreading 
the stresses over a much larger area. Consequently, the magnetic field 
broadens the interaction structure not by a large increase in the 
nonradial flow velocities (which are nearly the same in the HD and MHD 
cases) but by distributing these deflections over a greater range. 

For initially steep-sided streams. Figs. 2 and 3 demonstrate that 
dimensionality and magnetic effects act roughly on a par — neglect of 
either factor in a simplified model leads to gross misrepresentation of 
the evolutionary state Of the stream near earth orbit. We do note, 
however, that inclusion of the IMF declines in importance (relative to 
the dimensionality) for streams having broader initial transitions 
between fast and slow flows (e.g., were the input variations more 
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nearly sinusoidal) , In that case, the physical scale of the 
interaction region becomes so large that, even at the elevated 
magnetosonic speed, there is insufficient time for pressure signals to 
traverse the structure before it transits a substantial radial 
distance. | 

Thus the relatively economical 2-D MHD 'description emerges as the 
preferred vehicle for studying large-scale dynamics inside 1 AU and for 
experimenting with new physical processes (e.g., conduction, waves) in 
the structured solar wind. However, only the full 3-D formulation 
provides information on meridional flows, which are directly related to 
the unseen solar wind structure far from the ecliptic plane and which 
may have significant influence on the average angular momentum 
transport. 
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INVESTIGATION OF SECTOR BOUNDARY FINE STRUCTURE BETWEEN 0.3 AND 1 AU 

K.W.Behannon and F.M.Neubauer 
Inst, flir Geophysik und Meteorologie der Techn. Univ. 

3300 Braunschweig, F.R.G. 

Abstract Sector boundary regions as seen in magnetic field data taken 
by Hel ios 1 between 0.3 and 1.0 AU during its primary mission were gen- 
erally complex in structure, consisting of clusters of fine-scale di- 
rectional transitions. The total field magnitude typically decreased in 
association with these directional discontinuities. The results of min- 
imum variance analysis of 103 well-defined transitions show that the 
transition surface had a widely variable, at times perhaps oscillatory, 
local orientation. In 75% of the cases, the tilt of the surface relative 
to the ecliptic plane was found to be > 450 . The transitions were gen- 
erally most similar in magnetic characteristics to tangential discon- 
tinuities, although the results are not inconsistent with the observation 
of a small but nonzero normal component in at least half the cases stud- 
ied. 

1. INTRODUCTION 

In spite of a great deal of work on the macroscopic aspects of mag- 
netic sector structure over the past 14 years, very little has been 
done on the fine structure of the sector boundaries themselves. Studies 
to date (Smith, 1972 ; Bavassano et al . , 1976 ) have presented results 
based on a relatively small number of boundary observations, A survey 
over a large number of observed boundary transitions shows that such 
transitions are usually more than just a single simple reversal of mag- 
netic field polarity. They are often very complex in structure, con- 
sisting of a series of abrupt changes in field direction together with 
depressions in field magnitude. Thus, on the large scale, one is led 
to speak of "sector transition regions" in the observational data that 
consist of clusters of directional discontinuities and are also periods 
of relatively low and variable field strength. 

This report presents some of the preliminary results from a compre- 
hensive study of the sector transitions observed by the T.U, Braunschweig 
fluxgate magnetometer on Helios 1 (Musmann et al . , 1975). Specifically, 
the data are from the primary mission period, December 1974 to April 
1975, during which the spacecraft covered the heliocentric distance 
range 1.0 £ r j< 0.3 A. U. The results to be described are from the min- 
imum variance analysis of primarily 8 sec. average magnetic field vec- 
tors. 

Before analysis results are presented, examples of data on several 
time scales, including the high resolution data (4 vectors/ sec ) , will be 
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shown to illustrate both the large scale complexity and the fine-scale 
structure of sector transition regions. Quantitative results will then 
be described which address specifically 

(1) fine-scale boundary orientation, and 

(2) the type of discontinuity that best describes the fine-scale bound- 
ary. 

There were two major magnetic sectors during the Helios 1 primary 
mission period, with at least two very narrow sector regions also visi- 
ble (Neubauer, 1978). In the present study, eleven data intervals of 
from one to five days each have been analyzed. Within those intervals, 
a total of 103 fine-scale transitions were found which 

(1) appeared to be sector boundary-associ ated by visual inspection of 
the data; 

(2) had acceptably wel 1 -resol ved variance ellipsoids (each case satis- 
fied the cri, 2 rion where ^2 and are. the intermediate 

and minimum eigenvalues, respectively, determined in the minimum 
variance analysis); and 

(3) had angles of rotation in the discontinuity plane ai ^ 120° to dif- 
ferentiate them from ordinary directional discontinuities in the 
solar wind, since > 90% of such discontinuities are found to have 
to < 120° (Burlaga et al . , 1977 ; Lepping and Behannon, 1977 ). 

2. EXAMPLES OF TRANSITION REGION DATA 
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On the macroscale, the sector boundary which separated the two major 
magnetic sectors appeared to be a relatively simpl e-structured boundary. 
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complex in structure du- 
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and simple and sharply 
defined on the second and 
fourth recurrences. The 
complexity observed on 
pa.sses one and three con- 
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field, but also of peri- 
ods of intermediate 


270 

^ 1 ; j 1 ; 1 

• 

♦ 190 
90 

1 ! ' 1 ' 1 
r-v 1 1 1 

i! ^ it'* t lv/ J *’ 1. ’w/yl*' ■ i_l i 

t 1 1 



a 

« 0 

-45 

ipi rn r, 1 \ r r _ j 

. . . . . It. .i. . .*1 i. .1 1 . . 1 , L , 1 < i 1 1 1 i 1 1 1 I 1 1 1 1 1 i t 1 < 1 1 > 1 *1 1 ) 1 1 1 1 t i t 

1 1 

ait.- - j-'- ••-7 V - “ 

1 1 1 1 i l_A 1 1 i > 1 i 1 t-L 

Fig 

6 12 16 6 12 16 6 
IS JIEC 16 DEC 

1. Helios 15-min. magnetic field 

12 16 
17 DEC 

measure- 


ments of complex sector boundary near I AU . 
Data are field magnitude F, azimuth angle <j) in 
ecliptic plane (0° towa rd .,^.u;p-r) > and elevation 
angle e above and below ecliptic plane. 


OfilGIW/lL- PAGE 

POOR 


270 
<t) 180 


RMS 

UT 8:00 


10 20 

17 DEC. 1974 


(non-spiral) field orientation. 

Figure 1 shows 15-min, averages of field data taken during the first 
period at 0.98 All. During the two days prior to the major change from 
"away" (-) to "toward" (+) orientation^ the direction angles show that 

the field tended to al- 
■ ternate between the "away" 

F ' ■ spiral direction and a 

(y) 10 ^ n direction intermediate in 

5- VI - azimuth between the "away" 

- M I I I I I I I I I I I )-t-M I and "toward" directions. 

' A tendency for alternation 

4) 180 - .. ■ 

90 - ' - between positive and neg- 

-< I 1 1 I I I n t t M - M M - i - i I t t-t-i-i I I !■■(- ative 9 angles can also 

45 ^ , - ^ 

® - 45 I _-=r: seen. These data sug- 

-I I 1 1 M M I 1 I I I I I I II I 1 I I I I I I I I I gest that a series of 

RMS ’2' 

M il possibly filamentary 

UT 8:00 10 20 30 structures was observed 

17 DEC. 1974 in the sector transition 
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to field magnitude and direction, RMS 

deviation of magnitude F is also given. also observed to occur 

after the main, macroscale 
polarity reversal. In all, 
13 individual transitions 
that satisfied the above- 
listed selection criteria 
were found during the 
period shown, some of them 
corresponding to the 
"filament" boundaries 
seen in the 15-min. aver- 
ages. 
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Fig. 2. Helios 8-sec. averaged data showing 
fine-scale sector transition. In addition 
to field magnitude and direction, RMS 
deviation of magnitude F is also given. 
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Fig. 3. Detailed (250 msec) measurements of positive to negative 
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during the relatively sharp direction change. 

Figure 3 shows the most detailed data for this transition. One sees 
in this case a small but persistent normal component, ( <B 2 >= 1 '4y) • 

The continuous change in magnitude during the transition results in the 
roughly linear rather than circular appearance of the hodograph. 

3. ANALYSIS RESULTS 

1. A major subject of interest in this investigation concerned the orien- 
tations of the fine-scale transition boundaries studied. Minimum variance 
analysis yielded an estimate of the normal vector orientation in each 
case. Figures 4a and b show examples of these results for two different 
boundary transitions in terms of the e,(j) coordinates of the tip of the 
unit normal vector. For each example, arrows mark the azimuthal angle 
values of the directions perpendicular to the average observed field 
azimuth angle, <{)g, and to the spiral direction, calculated from the 
mean solar wind speed. The points plotted have been coded according to 
the value of the eigenvalue ratio, X 2 /^ 3 » Previously defined. 

In Figure 4a the results for the transitions comprising the boundary 
region shown in Figure 1 are given. Although the set of normal directions 
in this case are spread over ranges of 110° in e and 99° in a tendency 

is seen in the majority of cases for recurrence of an azimuthal orienta- 
tion centered roughly half way (within + 20°) between the radial direc- 
tion and the direction 
perpendicular to the 
spiral field direction. 
Some oscillation of the e 
angle is also seen , al- 
though in only three 
cases were negative val- 
ues obtained. 

Figure 4b shows re- 
sults from a transition 
region which consisted 
on the macroscale of a 
negati ve-to-posi ti ve 
polarity change, with 
‘*’n 11 fine-scale transi- 

Fig. 4. Orientations in elevation and azimuth tions satisfying the 
of collective transition surface normal vec- , . • • • r, 

tors for sector transition periods (a) near selection criteria. Here 

1 AD and (b) at 0.52 AU. For symbol definition the normal orientations 
see text. 
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were scattered about an azimuthal direction roughly perpendicular to the 
average spiral direction. There was also an oscillatory deflection of 
the normals often to large elevation angles, both positive and negative. 
Such data are consistent with a model of the boundary in which the sur- 
face is rippled or corrugated. 

A statistical study over all 103 fine-scale transitions being con- 


siderjpd in this investigation shows that <e^> <sf 0°, with < 45° in 

75% r/f the cases and = 32° + 21°, which means that the fine scale 


:° + 21 ° 

boundary surface was predominately oriented at an angle > 45° to the 
ecliptic plane. The wide variability seen in fine-scale boundary orien- 
tation makes it difficult, if not impossible, to infer an average orien- 

00 ^ I tat ion or "tilt" for the macroscale 

current sheet from such observations. 
2. The second major point of interest 
in the investigation was the identi- 
fication of the type of discontinuity 
most consistent with the fine-scale 
transition observations. Figure 5 
shows the distribution of normal 
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component magnitudes, normalized in 


each case by the average total field 
magnitude, i.e., |B 1/B . As can be 

mli^ituderorLctoJ'i^a^Juio^ 82Vof the cases had 

surface normal vectors relative to _< 0.3, a significantly larger fraction 

average total field in each case. ^5 found in general for ordinary 

directional discontinuities in the solar wind (see review by Neubauer 
and Barnstorf, 1978). In particular, few cases were found which were 
consistent with "1 arge- | B^ | /B^" rotational discontinuities. On the other 
hand, by conservative estimate approximately half the cases yielded nor- 
mal component magnitudes which were larger than one would expect for 
a truly zero normal, when possible measurement and computational errors 
are considered. 


4. SUMMAR Y 

A study of the sector transitions observed during the primary mission 
of Helios 1 shows that such transitions were generally "complex" rather 
than "simple" in structure on the fine-scale. Results suggest that the 
boundary surface may at times have had a locally irregular or "corru- 
gated" structure, since 

a. the transition regions were observed to be composed primarily of clu- 
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sters of di sconti nui ti es i 

b. the tilt cif the discontinuity surfaces was found to be > 45° on aver- 
age; and 

c. coherent patterns of oscillation of the transition surface normal 
direction were sometimes seen. 

It is concluded that the fine-scale boundary is more near! y ''TD-1 i ke" 

"RD-like" on average on the basis of 

a. the generally small magnitude found for the component normal to the 
boundary surfaces, and 

b. the observation usually of "dips" in total field magnitude in asso- 
ciation with the directional discontinuities. 

For at least half the cases studied the results were not inconsistent 

with the existence of a small, nonzero normal component as required, 

for example, for reconnection. 
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